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We develop a general method to calculate Zeeman splittings of electrons and holes in semicon- 
ductor nanostructures within the tight-binding framework. The calculation is carried out in the 
electron-hole picture and is extensible to the excitonic calculation by including the electron-hole 
Coulomb interaction. The method is suitable for the investigation of quantum shape effects and the 
anisotropy of the g-factors. Numerical results for CdSe and CdTe nanostructures are presented. 



I. INTRODUCTION 

Controllability of spins in semiconductor nanostruc- 
tures, or quantum dots, has become one of the impor- 
tant subject to be investigated in recent years due to 
the novel field of spintronicsi and quantum informa- 
tion processing^ Manipulation of the spin depends cru- 
cially on fundamental spin properties such as the effec- 
tive Lande g-factors^ On the other hand it is known that 
quantum confinement gives rise to novel electronic and 
spin properties and the shape control of nanostructures 
has become possible recentlyA^ It is thus imperative to 
understand what kind of quantum effects on magneto- 
optical properties can be induced by changing the size 
and the shape of nanostructures. 

The size and shape effects on g-factors in nanos- 
tructures have been investigated within the k • p 
framework 6,7,8 as well as the tight-binding frameworkAifi 
Within the kp framework the quantum confinement and 
shape effects are taken into account by imposing vari- 
ous boundary conditions. However, for the smaller size 
nanostructures the atomistic effects may become more 
relevant, and it is more difficulty to capture those effects 
with k • p method. For example, the mixing between 
heavy-hole band, light-hole band, and the spin-orbit split 
band become stronger as the size of the nanostructures 
goes down. The symmetry can be reduced, depending on 
the particular shape of the nanostructure. On the other 
hand it has been proven that the tight-binding approach 
can capture the atomic nature and surfaces effects. It 
does not need to assume the nanostructures possessing 
any symmetry, and has been used to calculate the g- 
factors of CdSe nanocrystals with reasonable successi*^ 

In this work we develop a general method which is 
suitable for the investigation of quantum shape effects 
on both electron and hole g-factors. We emphasize that 
we are working in the electron-hole picture instead of 
the conduction-valence band electron picture. This is 
one of the major difference in comparison with some 
of the preceding work using tight-binding methodAiS 
The electron-hole picture is more relevant to the real 
experimental and the method can be extended to cal- 
culate the exciton Zeeman splitting in a straightforward 
manner. The magnetic field is assumed to be in the 
weak magnetic regime, in which the cyclotron length 
Ib = \Jficje\B\ is larger than the size of the nanos- 



tructure. The quantum shape effects are investigated 
in detail. In the following we will first describe how to 
setup the Hamiltonian in clcctron-holc picture including 
the effect of vector potential associated with the external 
magnetic field. We then discuss in detail how to extract 
Zeeman splittings based on this Hamiltonian. Numer- 
ical results for CdSe and CdTe nanostructures will be 
presented and discussed. In the end we will discuss the 
effect of Coulomb interaction which is omitted in this 
work. 



II. METHOD 

The first step in the whole analysis is to obtain the 
single particle eigenfunctions and eigenenergies based on 
a semi-empirical tight-binding method, resulting in a 
Hartree-Fock ground state. The spin-orbit interaction 
is included in the tight-binding Hamiltonian. The dan- 
gling bonds are truncated. An electron-hole transfor- 
mation over the ground state is employed to obtain the 
Hamiltonians in the electron-hole pictured The second 
quantization field operator is written as 

*w = 2:«. + (*t ) ) + (1 » 

-ix(fc£0+2>(fc$). 

where a\ a is the annihilation operator of local orbitals 
with collective site-orbital index I and spin index a, and 
dfe c (»fe„) is the annihilation operator of the conduction 
(valence) band electrons. The field operator expanded in 
local basis gives rise to the typical single-particle tight- 
binding Hamiltonian 

Ht.b. = ^ U<j,i'(j'a\ a j/ai'a', (2) 

la,V a' 

where the summation is restricted to the on-site and 
nearest-neighbor pairs in this work. We use the 
well-known semi-empirical values for the tight-binding 
parameters^ By diagonalizing this Hamiltonian the 
eigenenergies Ek , the relation between a/ CT and Ofc, , as 
well the relation between <pi(x) and <f>ki{x) can be de- 
rived. To transform the Hamiltonian to the electron-hole 
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picture we define the hole operator to be hk v = a k . In 
this picture the free Hamiltonian without external mag- 
netic field becomes 

Ha = E E *A e a k. ~ E E Kh\h K . (3) 

Note that Ek c > 0, Ek v < 0, and we have dropped the 
overall energy constant of the Hartree-Fock ground state. 

The external magnetic field will introduce two addi- 
tional terms to the Hamiltonian, giving rise to Zeeman 
splittings. The first contribution comes from the Zeeman 
Hamiltonian: 

H Z eeman = J d? X^^tWdOVB (b ■ S) l/V^); 

<J(T ! 

(4) 

where go — 2.0 is the free electron g- factor, \i& is the 
Bohr magneton, S — \a, and B represents the external 
magnetic field. It can be shown that in the electron-hole 
picture the Zeeman Hamiltonian becomes 

Hzeeman = ^ 9°f Jl B e l c e k' e E E ' ^'/fe^'C 5 ) 

k c k' I uu' 

£ gofiBhl, h kv ]T E ft J ■ S aa , f k ,„, , 
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where we have defined f l k cr = J d?x<f)*(x)4>k iCT (x). Note 
that one should be careful about the indices of the hole 
operators. Within the tight-binding framework, sim- 
ply means the coefficients of the eigenfunctions in terms 
of local orbitals. The second contribution comes from the 
vector potential associated with the external magnetic 
field. It is known that the effect of vector potential A 
can be incorporated into the single-particle tight-binding 
Hamiltoniania by modifying the Hamiltonian to be 

H A = E t W>r' e ~ iXu ' a h,l> a l'<r'> ( 6 ) 

i(7,IV' 

where 

X U ,(B) = ~ C 1(B). dr. (7) 

Transforming the Hamiltonian into electron-hole picture, 
one finds that 

H A = H + E{B) (8) 

+ E e W E ^ i e ~ lXl " - !) ftJka- 

k c k' c laV <j' 
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where 

E & - E E ( e ~ lXl " - l ) fluL- ( 9 ) 
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represents a small overall energy shift of the ground state 
due to the vector potential. It should be noted that 
during the transformation we have dropped the terms 
which do not separately conserve the electron and hole 
numbers^ Note also that there is no magnetic field in- 
duced coupling between electrons and holes within this 
approximation. 

In this work the magnitude of the magnetic field is 
restricted to less than 15 Tesla. For an uniform field 
\B\ < I5T and size of nanostructure less than 100 A in 
average diameter, it is easy to estimate that X(B) is at 
the order of I0 -2 to 10 -3 , which means e x — I « X. 
As a result each term in Ha is roughly linear in B for 
the parameter range we are interested in. To summarize, 
the two-particle electron-hole Hamiltonian including the 
effect of external magnetic field is represented by 

H c h = -^Zeeman + -^A, (10) 

where Hzeeman is defined by Eq. [5] and Ha is defined by 
Eq. El 



III. ZEEMAN SPLITTINGS 

In this section we will briefly discuss how to extract 
Zeeman splittings based on Eq. ^3 starting from the 
single particle Hamiltonian Eq. [21 The first step is 
to find the single particle eigenenergies and eigenfunc- 
tions of Eq. |21 Once the single particle eigenenergies 
and eigenstates are known, the Hamiltonian in electron- 
hole picture, H c h = -ffzeeman + H\, can be readily con- 
structed. In this work we have used exact diagonalization 
for nanostructures with less than 777 atoms. It will be 
shown later that the full spectrum of the single particle 
Hamiltonian is not necessary if we are only interested in 
the Zeeman splitting of states within some energy win- 
dow. Hence for larger nanostructures one may use Lanc- 
zos or other methods to evaluate the single-particle states 
within some energy window, to reduce the computational 
demand. 

It is evident from the Hamiltonian, -ff e h, that the exter- 
nal magnetic field induces coupling between all zero-field 
eigenstates. In principle in order to evaluate the eigenen- 
ergies at non-zero magnetic field one has to diagonalizc 
the whole Hamiltonian, including all the electron levels or 
hole levels. However in nanostructures, quantum confine- 
ment effects give rise to a discrete density of states, lifting 
the degeneracies. In particular for the valence band, the 
heavy-hole light-hole degeneracy is usually lifted. As a 
result, typically the only degeneracy left is the Kramers' 
degeneracy. For a given a Kramers' doublet, the exter- 
nal magnetic field will remove the Kramers' degeneracy 
and introduce a Zeeman splitting for the doublet. For 
small external magnetic field the coupling to other dou- 
blets can be neglected compared to the coupling within 
the doublets, since there is a finite energy seperation be- 
tween different doublets. It is thus convenient to define 
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an intrinsic Zeeman splitting, which is determined by re- 
stricting the Hamiltonian to the Hilbert space spanned 
by one particular Kramers' doublet. 

As stated near the end of last section, for the magnetic 
field we are interested in, an intrinsic Zeeman splitting 
is nearly linear in B. We define the intrinsic g-factor as 
the ratio between the intrinsic Zeeman splitting and the 
magnitude of the magnetic field (up to a sign). When the 
intrinsic Zeeman splitting is smaller than the inter-level 
spacing the effects from nearby doublets are small, and 
the intrinsic Zeeman splitting is close to the true Zeeman 
splitting. However, when the intrinsic Zeeman splitting 
becomes comparable to the inter-level spacing the effects 
from nearby doublets become important. The intrinsic 
Zeeman splitting may substantially deviate from the true 
Zeeman splitting. In this case, it becomes necessary to 
include nearby doublets in order to find the true Zeeman 
splitting. 

We have numerically verified that typically less than 
5 nearby doublets are needed to be included. Inclu- 
sion of further states does not change the Zeeman split- 
ting. Since we are interested in the states nearby the 
band edge, we only need to find the single-particle lev- 
els near the band edge, eliminating the need to find the 
whole spectrum and eigenfunctions. When more than 
one Kramers' doublet is included in the calculation, the 
Zeeman splitting may become strongly non- linear in B, 
or level (anti-)crossing might occur. In those cases the 
g-factor becomes ill defined, but one can still define an 
intrinsic g-factor based on the intrinsic Zeeman splitting. 

The Zeeman splitting can be consider to be always pos- 
itive, while we usually associate a sign to the g-factor. ft 
is thus important to determine the sign of the g-factor in 
a consistent way. Due to the spin-orbit interaction and 
quantum confinement, there is a difference in spin content 
between the free carriers and the carriers in nanostruc- 
tures. In this work we use the intrinsic Zeeman splitting 
to determine the sign of the (intrinsic) g-factor. Basi- 
cally one has to find the correct zeroth order eigenstates 
within the Kramers' degenerate space, and use the spec- 
trum weight of zeroth order eigenstates to connect to the 
real spin. 

As an example to illustrate the general numerical pro- 
cedure to assign the sign of the g-factor, consider the case 
in which the external magnetic is pointing along in the 
z-direction. The first step is to find the eigenfunctions 
of -ff e h when a very small magnetic field is added. One 
then calculates the spectrum weight of the two eigenfunc- 
tions \r]), \n') and connects them to the real spin using 
their major spin component, for example, |ry) — > |s+) 
and | n 1 ) — > |s— ). Now if a zeroth order eigenfunction \n) 
moves toward higher energy as one increases the exter- 
nal magnetic field, the g-factor is assigned to be positive. 
In this analysis one must ensure that the spin quanti- 
zation axis aligned with the external magnetic field. As 
a result an unitary transformation on the eigenfunction 
might be needed when the magnetic field is not along the 
z-direction, the typical spin quantization axis. This as- 
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FIG. 1: Spectra weight of single particle eigenstates. 



signment may also become impractical when it becomes 
difficult to identify the major spin component of zeroth 
order eigenstates. In the results described below, we have 
found it difficult to identify the major spin component for 
some of the hole states. 



IV. RESULTS 

In the bulk limit, the conduction band is S = \ like 
and the valence band is J = § like. We calculate the 
spectra weight of zero-field single particle eigenstates in 
terms of local orbitals. In all calculations we adopt an 
coordinate system in which the z-axis is parallel to the 
c-axis of the wurtzite structure. We find that for the 
conduction band, around 90% of the spectra weight of 
the band edge state is in S = h, which increases as the 
size of nanostructures increases. However for the same 
nanostructure, the S — ^ component decreases as we 
move toward higher energy condunction electron states. 
For the valence band, we find that the band is J = § 
like within a fair large energy window near the valence 
band edge , containing around 90% of the spectra weight. 
However the mixing between heavy-hole and light-hole is 
strong and the mixing is sensative to the size and shape 
of the nanostructure. It is thus improper to assign the 
label of heavy-hole or light-hole to the holes states in 
those nanostructures. In Fig^we plot the spectra weight 
of zero-field eigenstates near the band edges for the 777 
atoms CdTe and CdSe nanostructures. Nanostructures 
containing different number of atoms show qualitatively 
similar behavior. The main features to be noticed are 
the reduction of S — ^ (S = f ) components for the elec- 
tron (hole) states as we move away from the conduction 
(valence) band edge, and the strong heavy-hole light-hole 
mixing. 

For zero field energy levels, the band edge states show 
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FIG. 2: Electron Zeeman splitting for 777 atoms CdSe and 
CdTe nanostructures. 
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FIG. 4: Non-linear Zeeman splittings, (a) Electrons in 561 
atoms CdSe nanostructure. (b) Holes in 576 atoms CdSe 
nanostructure. (c) Electrons in 561 atoms CdTe nanostruc- 
ture. (d) Holes in 450 atoms CdTe nanostructure. 
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FIG. 3: Hole Zeeman splitting for 777 atoms CdSe and CdTe 
nanostructures 
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FIG. 5: Intrinsic g-factor of first ten electron and hole states 
for CdSe nanostructures at various sizes. Filled circles repre- 
sent g z and empty circles represent g x . 



redshifting as the number of atoms in nanostructures in- 
creases. The valence band edge state occurs as heavy- 
hole like state for some of the nanostructures while oc- 
curs as light-hole like state for other nanostructures, even 
when there is only minor difference in the numbers of 
atoms contained in nanostructures. This result indicates 
that the shape of the nanostructure has profound impact 
on the spin content of the hole states. No degeneracy ex- 
cept for the Kramers' degeneracy is observed, which we 
attribute to the strong quantum confinement effect of the 
small nanostructures. In Fig[2l (FigEJ we plot the mag- 
netic field dependence of several electron (hole) energy 
levels for 777 atoms CdSe and CdTe nanostructures. The 
magnetic field is pointing along the z-direction, which co- 
incides with the c-axis of the wurtzite structure. As seen 



in these figures, the Zeeman splitting is nearly linear in 
B. Magnetic field dependence for most of the electron and 
hole levels calculated in this work show qualitatively the 
same behavior, while quantitatively the slope of Zeeman 
splitting varies, giving rise to different g-factors. How- 
ever some of the electron and hole levels show strongly 
non-linear dispersions. In Fig. 0] we plot the Zeeman 
splitting for some of those levels. The non-linear disper- 
sions occur when the intrinsic Zeeman spitting is close to 
the inter-level spacing, giving rise to level anti-crossing, 
or when the intrinsic Zeeman splitting is small compared 
to the splitting induced by inter-doublets couplings. 

In Fig[Sl(Fig||)J) we plot the intrinsic g-factor of first ten 
electron and hole states for CdSe (CdTe) nanostructures 
at various sizes. We found that g x rj g y ^ g z and g y is 
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FIG. 6: Intrinsic g-factor of first ten electron and hole states 
for CdTe nanostructures at various sizes. Filled circles repre- 
sent g z and empty circles represent g x . 



not plotted for clarity. In Table [I] we summarize the size 
and aspect ratio of the nanostructures used in FiglSJand 
FigED The aspect ratio is defined to be the ratio between 
L z and the average diameter ^/L x L y in x-y plane. The 
prominent features are the redshift (blueshift) of electron 
(hole) levels and anisotropy of g-factors. Quantitatively 
we found it difficult to identify a simple scaling relation 
between the g-factor and size or aspect ratio of nanos- 
tructures using the data shown here and other data for 
smaller size nanostructures. The size (L x , L y , L z ) and 
the aspect ratio are average quantities which are used to 
characterize the nanostructures. For higher dimensional 
systems, these average quantities have been successfully 
used to characterize the g-factors, since the detail of the 
boundary is less significant. However for small size nanos- 
tructures the detail of the structure is important in deter- 
mining the energy levels, wavefunctions, and hence the 
g-factors. For example, it is difficult to modify the as- 
pect ratio along without altering substantially the detail 
shape of the structure. We speculate that this is the main 
reason to make it difficult to identify a simple scaling re- 
lation between g-factor and other averaged quantities. It 
is thus advantageous to use tight-binding method, which 
captures the detail of the structure, instead of effective 
mass approximation type methods. 

Qualitatively some observations can be made. First 
since the g-factor is not a smooth function of the 
eigenenergy, nearby states might possess very different 
g-factors. The g-factors also change substantially as the 
size changes. We observe that some energy levels move 
toward higher energy while some levels move toward 
lower energy as one increases the aspect ratio. Similar 
dependence has been recently observed in the calcula- 
tion of electronic states of CdSe quantum rods^. We 
also observe that the electron g-factors are less sensitive 
than the hole g-factors to the aspect ratio. This can be 



TABLE I: Size and aspect ratio of the nanostructures. 



Number of atoms 


450 


561 


758 


768 


777 


\J L x L y (A) 


26.76 


22.85 


34.55 


27.42 


20.44 


L* (A) 


21.88 


35.88 


35.88 


39.38 


49.88 


Aspect Ratio 


0.82 


1.64 


1.04 


1.44 


2.44 



attributed to the fact that the S = \ like behavior of the 
electron is insensitive to the aspect ratio while for the 
J = i like holes, the mixing of heavy-hole and light-hole 
is sensitive to the aspect ratio. 



V. SUMMARY AND DISCUSSION 

In summary we have developed a general method to 
evaluate electron and hole Zeeman splittings in semi- 
conductor nanostructures within a tight-binding frame- 
work. The calculation is carried out within the electron- 
hole picture instead of the conduction- valence band elec- 
tron picture. Hence the scheme can be readily extended 
to the exciton Zccman splitting calculation by includ- 
ing the electron-hole Coulomb interaction. The results 
found here are qualitatively similar to the results in Rcf 
9, in which only band edge states are calculated. The 
quantitative difference may be attributed to the approx- 
imation used during the transformation to the electron- 
hole picture. In most cases we observe nearly linear Zee- 
man splittings induced by the external magnetic field. In 
some cases we also observe non-linear Zeeman splittings 
which are due to the interlevel coupling between nearby 
Kramers' doublets. Those results are qualitatively simi- 
lar to what are obtained in RefQ. 

We find that the behavior of electron Zeeman splittings 
of nanostructures with different aspect ratio does not 
show drastic difference. This is in contrast to the dras- 
tic difference (change of sign) between spherical quan- 
tum dot (SQD) and semi-spherical quantum dot(SSQD) 
found in RefQ- It should be noted that in RefQ SSQD 
is treated by imposing boundary condition on envelope 
function, requiring envelope function to be vanishing 
in the upper hemisphere. This boundary condition is 
very different from the boundary condition used here for 
nanostructures with different aspect ratio. In other word, 
the idea of SSQD in Ref □ is different from the idea of 
quasi-spherical regime for the elongated nanostructure 
with varying aspect ratio in other referencesA^iifl 

There may be strong difference between magneto- 
optical properties of free electron-hole pairs and ex- 
citons, in which the Coulomb interaction is included. 
Two key questions to be addressed are the effect of 
the Coulomb interaction and how to experimentally ob- 
serve electron, hole, and exciton Zeeman splittings us- 
ing magneto-absorption spectroscopy. To answer the first 
question one has to estimate how strongly the electron- 
hole Coulomb interaction modifies the free electron-hole 
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pair picture. In Refill, a restricted basis, single exci- 
tation configuration interaction representation method 
is used to calculate the exciton fine structure for CdSe 
nanostructures. It was shown that about 18 electron and 
24 hole single-particle states needed to be included to get 
a reasonable accurate exciton energy. We expect that 
roughly the same number of states will be needed in or- 
der to accurately evaluate the exciton Zeeman splittings. 
Since our results indicate that the g-factor of electrons 
and holes may vary substantially from one state to an- 
other, we expect that the exciton g-factor may be very 
different from the simple sum of the g-factor of free elec- 
trons and holes. Starting from the general framework 
developed here, the effect of Coulomb interaction can 
be included via configuration interaction methoo"2*A& or 
time-dependent tight-binding method^ 

To shed some light on the second question we have to 
consider the problem of optical orientation in nanostruc- 
tures. In bulk or quantum well materials, one typically 
assumes that the conduction band edge is s-like while 
the valence band edge is heavy-hole like. In this case 



the optically active excitons only consist of heavy-holes. 
Due to the angular momentum conservation, the opti- 
cally excited electron and heavy-hole spin will be aligned 
in particular configuration depending on the polariza- 
tion of the optical pulse. However, in nanostructures the 
strong heavy-hole light-hole mixing and the Coulomb in- 
teraction break down this naive picture, as it becomes 
difficult to determine the spin configuration of the con- 
stituent electrons and holes of the excitons. A realistic 
calculation of the optical absorption oscillator strength 
might be needed in order to fully understand what is 
actually measured in the experiment. 
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